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ABSTRACT 

Impulsive solar energetic particles (SEP) bursts are frequently observed in 
association with so-called eruptive flares consisting of a coronal mass ejection 
(CME) and a flare. These highly prompt SEPs are believed to be accelerated 
by the flare rather than by a CME shock, but in the standard flare model the 
accelerated particles should remain trapped in the corona or in the ejected plas- 
moid. In this case, however, the particles would reach the Earth only after a 
delay of many hours to a few days. We present a new model that can account 
for the prompt injection of energetic particles onto open interplanetary magnetic 
flux tubes. The basic idea underlying the model is that magnetic reconnection 
between the ejection and external open field allows for the release of the ener- 
getic particles. We demonstrate the model using 2.5D MHD simulations of a 
CME/flare event. The model system consists of a multipolar fleld with a coro- 
nal null point and with photospheric shear imposed at a polarity inversion line, 
as in the standard breakout model, but we include an isothermal wind in this 
system, which leads to the presence of open flux neighboring the eruption. We 
present the results of the simulation and compare these results with the standard 
observations of impulsive SEPs. We also discuss the implications of the model 
for observations and for further calculations. 
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Introduction 



Eruptive flares are well known to be the drivers of the most destructive space weather at 
Earth and in the heliosphere. Among the more hazardous forms of space weather, especially 
for human space flight, are the intense Solar Energetic Particle (SEP) bursts associated with 
fast coronal mass ejections (CME) /eruptive flares. SEPs are also important as a basic physics 
phenomenon, because particle acceleration is observed to occur throughout astrophysical 
and solar system plasmas. Two candidate mechanisms for particle acceleration in eruptive 
flares have been proposed: shock acceleration by the shock wave driven by a fast CME, 
and some type of Fermi, stochastic or electric field acceleration produced by the magnetic 
reconnection that drives the flare. These two mechanisms are believed to account for the 
classic observation that SEP events ap pear to be of two types: gradual events due to the 
CME-shock accele ration ( lReameslll999l ). and impulsive events due to the flare reconnection 
(ICane et al.lll986[ ) (hereafter called flare- accelerated particles). 



In order to reach the Earth, energetic particles must be injected from the acceleration 
site to interplanetary magnetic field (IMF) lines connected to the Earth, and along which they 
can propagate. For energetic particles accelerated by a shock ahead of a CME, the injection 
occurs when energetic particles reach velocity high enough to escape from the shocking 
region. The escaping particles are therefore injected from the shocked regions directly onto 
the open IMF lines. The particle injection from the shock to the interplanetary medium has 
been confirmed by a m ulti-instrument analys is of a particular SEP event detected at different 
longitudinal positions (IRouillard et al.ll201ll ). Also, this injection process at the CME-driven 
shock can theoretically explain the structure of the time profile of energetic particle fluxes 
( lRodrfguez-Gasenll201ll ). 



In contrast, the interplanetary injection of flare-accelerated particles that produce im- 
pulsive SEPs is far from straightforward. Flare acceler ation is be lieved to occur low in the 
corona, in the closed magnetic field of active regions (ILinl |2005| ). implying that energetic 
particles do not have direct access to open interplanetary magnetic field. In order for parti- 
cles to escape onto open flux tubes, the magnetic reconnection that accelerates the particles 
should involve both closed and open magnetic flux. This reconnection between the closed 



corona and the open IMF is referred to as "interchange reconnection" (jPariat et al.l 12009 



Edmondson et al 



20091 : iMasson et al.ll2012l ). Interchange reconnection typically occurs at 



null points present in magnetic configurations showing a transition between open and closed 
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magnetic field f lTitov et al.ll201l[): for example, the structures that are identified as hel met 
streamers, pseudo-streamers ( Wang et al. 2007al lbl) and coronal jets f Cirtain et al. 2007 ). 



Although particles can be accelerated at an intercha nge reconnection site and therefore 
have direct access to the open interplanetary medium (IMasson et al.ll2012l ). most of the 
impulsive SEPs are associ ated with eruptive f lares, which are usually not expecte d to involve 
inter c hange reconnection fjYashiro et al.ll2004r). The standard 3D CSHKP model (jCarmichael 



19641 : ISturrocklll966l : lHirayama]ll974l : iKopp fc Pneumanlll976l ) describes the global evolution 



of a solar eruption consisting of a flare and a CME. Figure [T] illustrates the two main phases 
of the CSHKP model and its implications for fiare-accelerated particles. In a dipolar active 
region, the two polarities are connected by a magnetic arcade that overlies an initially stable 
magnetic flux rope, corresponding observationally to a sigmoid/filament /prominence. When 
the flux rope becomes unstable, it rises in the corona and the overlying field lines stretch 
outward. A cu rrent sheet forrns within the arcade belo\y the flux rope and above the polarity 
inversion line (ILin et al.ll2005l l2008l : iReeves et al.l 120081 ). Eventually, magnetic reconnection 
begins in this current sheet , st art ing at low altitudes, then successively at higher altitudes 
(see review by iForbes et al.ll2006l ). The resulting reconnected field lines wrap around the 
flux rope, building it up further, and also close down to form the post-flare loops. Flare- 
accelerated particles are expected to be energized at the reconnectio n site below the flux 
rope, and thereafter to propagate along the reconnected field lines ([Gorbachev fc Somov 
19891 ). Hence, the flare-accelerated particles should be injected into the post-flare loops and 
the twisted field lines of the flux rope. In a 3D geometry, both footpoints of the flux-rope 
field lines are anchored to the solar surface, implying that energetic particles injected in the 
reconnected twisted field lines are trapped in the flux rope and do not have access to the 
interplanetary open magnetic field. 

The problem with this picture is that impulsive SEPs are observed at Earth a few days 
before the CME and ha.ye been clearly identified as accelerated during the imp ulsive phase 
of the flare ( Klein et al.lll999l : iMiroshnichenko et al.ll2005l : iMasson et al.ll2009al ). Therefore, 
the flare-accelerated particles somehow find a way to escape to the open interplanetary mag- 
netic field. In analogy to the interchange reconnection in a classical null-point topology, 
the coupling of the closed CME field to the open interplanetary field through magnetic re- 
connection may po ssibly provide a path for the escape of the initially trapped particles. 
Cohen et al.l (120101 ) suggested that during CME propagation, the CME flux rope interacts 
with the ambient magnetic field, affecting significantly the propagation of the CME. Multi- 
wavelength (white light, soft X-ray, and radio) analyses of some specific events showed that 
energetic electrons are injected along open magnetic field from the edge of the CME, sug- 
gesti ng that magnetic reconnection occurs between the CME and the open ambient magnetic 
field jMaia et al.ll2007t iHuang et allboili bemoulin et al.ll2012h . 
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It seems evident from Figure [T] that even with the presence of nearby open flux, an 
eruptive event in a purely bipolar region is highly unlikely to produce impulsive SEPs. On 
the other hand, we have argued that a multi-polar magnetic topology is essential for a fast 
CME/eruptive flare fjAntiochoslll998l ) and many observations have shown t hat eruptive flares 
are g enerally associated with complex, multi-polar active regions (e. g. iPatty fc Hagyard 
19861). The key featu re of eruption in a multi-polar topology, as in the breakout model 
( jAntiochos et al 1999 ). is that reconnection occurs between the fleld overlying the erupting 
flux and neighboring flux systems. Note, however, that this reconnection is physically distinct 
and spatially separated from the flare reconnection that accelerates the impulsive SEPs; 
consequently, the question remains as to whether a breakout eruption with nearby open flux 
can lead to the efficient escape of flare-accelerated particles. 

To answer this question we performed a 2.5D MHD simulation of an axisymmetric 
quadrupolar magnetic configuration, opened by the solar wind and forced by photospheric 



shearing motions, which trigger an eruption through the breakout model (jAntiochos et al. 



19991 ). We investigate the detailed dynamics of the magnetic field during the formation of 
the CME and its propagation into the open interplanetary medium. Previous related nu- 
merical investigations have obtained fast eruptions from a comp letely closed, static corona 



f lAntiochos et al.lll999l : iMacNeice et al.ll2004 



Karpen et al.ll2012[ ) or have yie lded only slow, 



Zuccarello et al. 


2008, 


2009. 


2012; 


Soenen et al. 


2009 


)• 



ported a fast breakout eruption from below a helmet str eame r, but only when n ear-sonic 



footpoint motions were imposed. Both ICohen et al.l (120101 ) and iLugaz et al.l (1201 if ) assumed 



a strongly out-of-equilibrium fiux rope, inserted into the initial state below a helmet streamer 
opened by the solar wind, to drive a high-speed eruption. In contrast, we obtain a fast CME 
with vigorous reconnection occurring between the ejected flux rope and an ambient open 
fleld, in response to gentle driving by subsonic footpoint motions at the base of our gradu- 
ally evolving conflguration. 

After describing the numerical model and initial conditions in ^ we present our results 
on the magnetic-reconnection dynamics during CME initiation (§3] and §1]), eruption (^, 
and propagation into the solar wind In §71 we describe the topological evolution 

resulting from the complex dynamics of the CME interacting with the solar wind. Finally, the 
implications of our results for energetic-particle injection into the heliosphere are discussed 
in m 
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2. Model description 
2.1. Equations, grids, and boundary conditions 

The simulations we re performed usinp; the Ad aptively Refined Magnetohydrodynamics 



Solver (ARMS; see e.g. iDeVore fc Antiochosll2008l ). solving the following ideal MHD equa- 



tions in spherical coordinates: 

|^ + V-(pu) = 0, (1) 

^ + V-(puu) = (l/47r)(VxB)xB-VP + pg, (2) 

— -Vx(uxB) = 0, (3) 

where p is the mass density, u the plasma velocity, B the magnetic field and P the pressure, 
and g = GM^r/r'^ the solar gravitational acceleration. We assume a fully ionized hydrogen 
gas, so that the plasma pressure P = 2{p/mp)kBT, where T is the temperature. For sim- 
plicity, we assume also that the temperature is constant and uniform, T = Tq, which allows 
high-lying field lines to be opened to the heliosphere by the solar wind while low-lying field 
lines can remain closed to the Sun, as occurs ubiquitously in the corona. Our objective in this 
paper is to simulate with high fidelity the changing connectivity of the near-Sun magnetic 
field, not to predict the detailed thermodynamic properties of the far heliospheric plasma. 
Thus, we adopted the simplest possible model that yields a solar wind and a dynamically, 
self-consistently determined boundary between open and closed coronal magnetic structures. 

The n umerical scheme is a finite- volume multi- dimensional Flux Corrected Transport 



algorithm (iDeVord Il99ll ). It ensures that the divergence of the magnetic field remains 
of the order of the machine accuracy, and that unphysical oscillations in all variables are 
prevented while introducing minimal residual numerical diffusion. The equations are solved 
using a second-order predictor-corrector in time and a fourth-order integrator in space. In 
conjunction with the fiux limiter, this accuracy is sufficient to inhibit numerical reconnection 
until any developing current sheets (e.g., at the pre-existing null point; §2.3p are compressed 
strongly to the grid scale by the ideally driven plasma fiows. The numerical resistivity then 
switches on to a small but finite value determined by the local fiow speed, grid spacing, and 
magnetic-field profile. 

The spherical computational domain covers the volume r G [l-R©, 125i?0], 6 G [— 7r/2, 7r/2] 
and (f) G [— 7r,7r]. The grid is stretched exponentially in radius r and spac ed uniformly in 



latitu de 6. ARMS uses the parallel adaptive meshing toolkit PARAMESH flMacNeice et al. 



2000l ) to tailor the grid to the evolving solution. In this simulation, the mesh was adaptively 
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refined and coarsened over six levels of grids, determined by tlie ratio of the local scale of 
the three components of electric-current density to the local grid spacing. Thus, the intense 
current sheets developing in the system are always resolved as finely as necessary on the mesh 
until the limiting resolution is reached. In addition, we imposed a maximum-resolution layer 
of grid cells at the inner radial boundary in order to resolve throughout the simulation the 
flows and gradients generated by the boundary forcing ( §2.41) . 

Both the inner and outer radial boundaries of our simulation are open to flows of mass, 
momentum, energy, and magnetic flux. We imposed line-tied conditions at the first radial 
cell inside the inner boundary, so that the magnetic field does not move tangentially except 
where forced by a prescribed flow velocity ( §2.41) . Over the three radial guard cells below 
the inner boundary, we kept the mass density p and pressure P fixed at their initial values 
( §2.2p . set the velocity v to zero, and applied a zero-gradient extrapolation to the three 
components of the magnetic field B from the interior values just above that boundary. At 
the outer radial boundary, we extrapolated both v and B using zero-gradient conditions, 
and applied a fractional multiplier to p and P that was set by the initial atmosphere and by 
the local radial grid spacing. The resulting inflow of material at the inner boundary evolves 
self-consistently in response to the changing solution in the interior and to the solar-wind 
outflow of material at the outer boundary. The polar boundaries of our domain were closed 
to all flows. We applied zero-gradient conditions there to p, P, f^, and S^, and reflecting 
conditions to f^, w^, Bq, and B^. 



2.2. Initial atmosphere 

In the steady state, the boundary between open and closed magnetic domains ultimately 
is determined by the force balance between the outward solar-wind pressure and the inward 
magnetic-field tension. To find this steady state for our starting magnetic configuration, we 
initialized the plasma using the spherically symmetric, isothermal, trans-sonic solar wind of 



Parker! (jl958[ ). The steady solution for the radial velocity v{r) of the atmospheric plasma 



flow at uniform temperature Tq can be expressed as 



^-p(i-^j = ;l-p(4-4-). (4) 

where Cg is the isothermal sound speed (c^ = IkBT^Im^ and r., = GM^mpj AksT^ is the 
sonic point. We assume a constant temperature To = 2 x 10® K, for which v = Cg = 
180 km s^^ at r = = 2.9Rq. This yields an acceptable solar-wind speed of 420 km s~^ 
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at 12O-R0. The inner-boundary mass density is a free parameter that we set to p{Rq) = 
1.63 X 10~^^ g cm~^, which yields reahstic /3 values (ratio of plasma thermal pressure to 
magnetic pressure) throughout the computational domain (see §2.31 and Figure [2]). 

The velocity v{r) computed from Parker's isothermal solution (Eq. |4]), together with the 
associated mass density implied by the steady mass- flux condition pvr"^ = constant, defines 
the initial state of the plasma throughout the numerical domain. We superimpose on that 
solution potential-field magnetic dipoles to create a quadrupolar magnetic configuration 
( §2.3p . Of course, this combined system initially is out of equilibrium, and must first be 
allowed to relax to a new quasi-steady state. We found that this required a relaxation 
time of about 2.1 x 10^ s, after which the kinetic and magnetic energies were essentially 
constant and no significant further evolution of the heliospheric plasma and magnetic-field 
distributions was observed. The left panel of Figure [2] shows a radial cut, starting at the 
polarity inversion line of the embedded active-region dipole, of the radial plasma velocity. 
The velocity plateaus at a speed of 500 km s"^ at a radius of 20Rq. 



2.3. Null-point topology 



Our initial magnetic configuration emulates the magnetic geometry of many solar erup- 
tions, here simplified to consist of a single dipolar active region embedded in one of the solar 
hemispheres. Such quadrupolar configurations have been shown to s uccessfully trigger fast 



CME eruptions through breakout initiatio n in 2.5D MHD simulations (lAntiochos et al. 



MacNeice et al.ll2004j : iKarpen et al.ll2012l ) and also in fully 3D situations ( iLynch et al 



1999 : 



2008 



2OO9I ). This study addresses the dynamics of a flux rope erupting near a coronal hole opened 



by the outward solar-wind flow acting on the background solar magnetic field. A global-scale 
dipole located at the center of the Sun, with a peak strength of 10 G at the photospheric 
poles, defines this background field. A smaller-scale magnetic dipole, representing the active 
region, is located in the northern hemisphere near the polar coronal hole. 

To simplify the problem still further for this initial investigation, we required the mag- 
netic field and the associated MHD solution to be invariant along the direction. We 
constructed a suitable 2.5D active-region magnetic field by assuming an azimuthally sym- 
metric toroidal ring of dipoles encircling the Sun. The resulting analytical expression for the 
vector potential of the dipole ring in Sun-centered cylindrical coordinates (^, 0, z) is (see the 
Appendix for details) 
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X 



AJ^ z) = 2Bo<PRq ^ 

U + ay + {z- hyf' [(e - ay + {z- by] 

[[{i + a)fi, + {z- E{k) - 2 [a/i, + {z - {K{k) + A;"^ - i^'lA;)]}} (5) 



where a and 6 are the ^ and z coordinates of the ring dipole; /i^ and /i^ are the ^ and 2; 
direction cosines of its magnetic moment; Bq is its field strength and d its depth below the 
surface; and K{k) and E{k) are complete elliptic integrals of the parameter 



{^ + ay + {z-hy 

To obtain the magnetic field created by the dipole ring in the domain volume, we then 
numerically differentiated the vector potential, B(r) = V x A(r). 

The location, depth and orientation of the dipole ring were chosen to yield a quadrupolar 
configuration haying a coronal null point w ith inner and outer spine lines and a fan surface 



( lAntiochos et al.l l2002l : iTorok et al.l |2009[ ). In the solar context, the outer spine can be 
either open to the interplanetary medium or closed to the photosphere, depending upon 
whether the active region is, respectively, embedded in a coronal hole or confined below 
large-scale magnetic loops. For this study, we adjusted the parameters to obtain an initially 
closed outer spine confined below a large-scale helmet streamer, after the system relaxed 
to its steady state. The left and middle panels of Figure [3] show respectively the initial 
field, at time t = s, and the magnetic configuration following the relaxation phase, at 
t = 2.1 X 10^ s. This initial null-point topology with a closed outer spine is characterized by 
its field strength Bq = 40 G, depth d = O.IRq below the photosphere, position at latitude 
27° from the equator in the northern hemisphere, and orientation (/ig, /i^) at a 45° angle from 
the southward tangent to the photosphere, pointing into the Sun. Magnetic separatrices pass 
through the null point and delimit the connectivity domains. These topological objects are 
labeled in the right panel of Figure Ei The dome-like fan separatrix surface bounds the inner 
connectivity domain, enclosed below the fan, and a closed outer connectivity domain, above 
the null point but within the helmet streamer. The fan surface is intersected at the null 
point by a singular spine separatrix field line: the inner spine is rooted below the fan, while 
the outer spine belongs to the closed outer connectivity domain. 

These dipole-ring parameters also yield a realistically variable plasma beta throughout 
the numerical domain. The right panel of Figure [2] shows a radial cut of /3 from the photo- 
sphere to 2Oi?0, starting at the polarity inversion line of the embedded active-region dipole. 



- 9 - 



One notices that /3 ^ 10~^ at the lower, photospheric boundary and /3 ^ 1 throughout 
the corona r < ARq, except near the coronal null point where the beta sharply peaks to 
a large value (> 10^). For r > 4i?Q, in the interplanetary medium, /3 rises gradually and 
continuously through the rest of the heliosphere out to the outer domain boundary, even- 
tually exceeding but remaining on the order of unity. Thus, our simulation is performed 
in a solar/heliospheric regime that is strongly field-dominated in the inner corona, but is 
plasma-dominated both at the null point and far out in the heliosphere. This contrasts with 



previous 2.5D simulations o f breakout configurations \vith a solar wind (Ivan der Hoist et al. 



20071 : IZuccarello et al.ll2008l . l2009l . l2012l : ISoenen et al.ll2009f ). in which the background and 
active-region field strengths were well below 10 G and the minimum plasma beta was 10~^ 
or larger. 



2.4. Boundary forcing 

The system is driven by photospheric motions applied to the magnetic flux confined 
below the fan's south lobe. Two photospheric flows are prescribed to displace field hues 
along the (j) coordinate in opposite directions on either side of the polarity inversion line 
(PIL) of the active region, thus forming a sheared arcade below the coronal null point. Each 
flow is defined by a cosine profile that avoids the formation of strong shearing gradients at 
the boundaries of the flows. Since the magnetic flux distribution below the south part of 
the fan is not symmetric, the photospheric flows on the two sides of the PIL do not extend 
symmetrically in the latitudinal direction. The velocity fields are applied over the ranges 
e G [24.3°, 28.1°] and 9 E [19.1°, 22.4°], respectively north and south of the PIL (right panel 
of Figure [3]). This photospheric forcing starts at t = 2.1 x 10^ s, after the relaxation phase. 
Hereafter, we consider the starting time of this forcing as the initial time of the simulation, 
t' = 0. The velocity fields are applied gradually by multiplying the spatial velocity field 
with a time- dependent cosine ramping function. The speed of the prescribed flows on both 
sides of the inversion line reaches a constant value after a time interval tpj^^^ = 1 x 10^ s, 
with a maximum driving velocity = 20 km s~^. This speed is greater than that of 

observed photospheric motions, for numerical convenience and efficiency. Nevertheless, the 
flows are still subsonic and highly sub-Alfvenic, with acoustic and Alfven Mach numbers of 
or der 10~^ and 10~^, r e spect ively. We note that these values are close to those employed 



by Ivan der Hoist et al.l (120091 ) in a 3D simulation of a slow streamer-blowout eruption; to 
obtain a fast CME, however, they had to impose footpoint flows that were faster by a factor 
of four. 



The imposed photospheric motions inject free magnetic energy and magnetic stress 
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into a force-balanced system. Electric currents are expect ed to deve l op along the variou s 
topological features, such as the null points and sep aratrices (lLowlll987l : lAlylll990l : lLaij|l993l ) . 
and within the sheared arcade flForbes et al.ll2006l and references therein). These structures 
are described below. 



Formation of volume and sheet currents 



The left column of Figure H] displays a 2D cut in the (r, 6) plane of the (f) component 
of the current density at times t' = 2.32, 2.58, and 2.74 x 10^ s after the shearing motion 
has reached constant velocity. Black and white correspond to the maximum values in the 
positive and negative directions, respectively. The white structure shows the current 
density developing within the sheared arcade. As the simulation evolves, the shearing motion 
induces the elongation and strengthening of the current structure, which occurred iri earlier 
simulations of the 2.5D breakout model (IMacNeice et al.l l2004l : iKarpen et al.ll2012l ). This 
distributed volume current is associated with the growing inward tension force exerted by 
the poloidal magnetic field {Br, Be), which balances the outward pressure force exerted by 
the shear-induced toroidal magnetic field component B^. 

T he applied shear ing motion also leads to the formation of a current sheet at the null 
point jAntiochosI [l996h . which IS the thin, quasi-horizontal black structure well above the 
white volume current within the sheared arcade (left column of Figure Hj). The right col- 
umn of Figure H] shows the evolution of the separatrices (dark blue field lines) and other 
magnetic-field lines involved in the dynamics of the simulation (other colors, see below in 
§1]). The shearing of the arcade causes its fiel d lines to bulge ou tward, compressing and 
deforming the separatrices and the null point (jPontin et al.l 120071 ). This leads to the mis- 
ali gnment of the inner and out e r spines as displayed by the dark blue field lines of Figure 
H JCalsgaard fc Nordlundlll994l : lAntiochos et al.ll2002l : iMasson et a h l2009bh. This she aring 



of the spines elongate s and intensifies the null-point current sheet flRickard fc Titovl Il996 



Galsgaard et al.ll2003l ). Magnetic reconnection is therefore expected to ensue first at the 



distorted null point. 



4. Reconnection at the null point 

The right column of Figure H] (see also the movie available in the online version) presents 
the connectivity evolution of selected magnetic field lines, which are plotted from fixed r and 
6 footpoint locations at each time. Gray lines represent magnetic field whose connectivity 
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does not change during the simulation, while the other colors represent specific connectivity 
domains. The red and orange field lines are initially closed below the fan's south lobe; 
their fixed footpoints are anchored in the negative (southern) polarity of the dipole ring. 
The yellow, pink, and green field lines are rooted in the negative polarity of the northern 
hemisphere and, respectively, are connected to the southern positive hemisphere, initially are 
opened to the interplanetary medium, and define the open fiux belonging to the northern 
coronal hole throughout the evolution. 

As the simulation evolves, magnetic reconnection takes place at the null point and the 
connectivity of the orange and yellow field lines changes. At t' = 2.32 x 10'^ s, the outermost 
orange line reconnects at the null point with a yellow helmet-streamer field line, which 
initially confines the closed outer spine. Subsequently, additional orange lines reconnect, 
transferring part of the fiux overlying the sheared arcade (represented by the red lines) 
toward closed magnetic domains located on both sides of the fan's south lobe (panels d & e 
of Figure H]). The yellow lines jump from the outer connectivity domain to the northern-most 
inner domain and form new loops below the fan's north lobe, while the orange lines pile up 
above the arcade connecting the two hemispheres, under the outer spine. 

In the right column of Figure HI the red and orange field lines belong to the same 
connectivity domain, but they reconnect with two distinct magnetic fiux systems. After the 
orange lines have entirely reconnected with yellow lines at the null point, the red lines start 
to reconnect there with the pink lines. Thus, the red field lines jump from the central inner 
connectivity domain to the outer connectivity domain. These sheared-arcade field lines, 
initially closed below the fan, now open into the heliosphere; meanwhile, the initially open 
pink lines close down, below the fan's north lobe (panel f of Figure H]). This exchange of 
connectivity between open and closed field lines corresponds to the interchange reconnection 
mode. 

During the transition from the classical closed/closed (orange/yellow) null-point re- 
connection to the interchange closed/open (red/pink) reconnection, the magnetic topology 
changes. During the closed/closed reconnection, the outer spine remains closed below the 
helmet streamer (panel e of Figure H]), while during the closed/open interchange reconnection 
the outer spine has been opened to interplanetary space (panel f of Figure H]). When the 
entire closed magnetic fiux (yellow) initially confining the null point has reconnected and 
been transferred to the closed domain of the fan's north lobe, the outer spine separatrix 
field fine is positioned at the boundary between the closed (red fiux) and open (pink fiux) 
connectivity domains. The outer spine is now opened in the interplanetary medium. This 
topological change implies that the null-point topology is now embedded in the northern 
coronal hole. Thus, the sheared arcade now is positioned to erupt into an open-field region. 



-VI- 



rather than into a closed-connectivity domain as usually studied. 

This interchange reconnection plays two crucial roles. First, it contributes to the flux 
removal, participating in the loss of equilibrium of the flux rope that is required in the 
breakout model of eruptions. Second, when the flux rope does erupt, its closed coronal flux 
is driven to interact and reconnect with the open interplanetary magnetic field. 



5. Flare reconnection and flux-rope eruption 

As the null-point reconnection progresses, the volume current within the sheared arcade 
becomes more distended, thinner, and more intense, due to the ongoing photospheric shearing 
motion (^. This evolution is evident in Figure [Hand accelerates strongly at the later times 
shown in Figure [51 When the volume current becomes very distended, its central section 
thins to a quasi- vertical current sheet that gets squeezed to the grid scale. The sheared 



arcade then reconnects in ternally, signaling the onset of the fiare-reconnection phase (IForbes 



2OOOI : iKarpen et al.ll2012l ). As in the schematic CSHKP model for eruptive flares (see ||T]and 
Figured]), two new flux systems are born: the post-flare loops, forming an arcade that closes 
down below the flare-reconnection site; and twisted field lines high in the corona, forming the 
ejected flux rope. Due to the azimuthal symmetry of the 2.5D simulation, the new flux-rope 
field lines are disconnected completely from the photosphere. 

The colored field lines plotted in Figure [5] are identical to those in Figure |H The field 
lines are over-plotted on a 2D cut in the plane (r, ff) of the component of the current 
density, j^, shaded in grayscale. The evolution of shows the evolution of the flux rope, 
while the colored field lines show the dynamics of the magnetic flux systems involved in the 
eruption. The six panels of Figure [5] display different times of the simulation - t' = 3.11, 
3.16, 3.17, 3.18, 3.20, and 3.26 x 10^ s, respectively - starting after the opening of the outer 
spine. (A movie is available in the electronic version of this article.) 

At t' = 3.11 X 10^ s, the structures are similar to those observed during the null- 
point reconnection (see Figure H]), except that the apex of the volume-current structure 
has acquired an ellipsoidal shape. The imminent formation of the disconnected flux rope 
suggests that the flare reconnection may have started already at this time. Examining our 
MHD simulation data at a high temporal cadence of 100 s, we looked for a magnetic 0-point 
in the numerical domain that would constitute the central axis of the flux rope. At time 
t' = 3.07 X 10^ s, an 0-point appeared at the top of the flare current sheet, confirming that 
the internal reconnection of the sheared arcade indeed had begun and that the flux rope was 
beginning to form. 
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Between t' = 3.11 x lO'' s (Figure [5ji) and t' = 3.18 x 10'' s (Figure ISH), the ellipsoidal 
current structure enlarges as the flare current sheet thins. This implies that sheared-arcade 
field lines are reconnecting at the fiare current sheet. The newly reconnected field lines wrap 
around those that reconnected previously, increasing the fiux-rope magnetic fiux. Meanwhile, 
the fiux transfer enabled by the null-point and interchange reconnections (§!]) induces a 
decrease in the downward magnetic tension of the overlying magnetic field . These changes 



in the force balance drive t he rise and eventua l take-off of the fiux rope (lAntiochos et al. 



19991 : iMacNeice et al.ll2004l : iKarpen et all 120121 ). 



In Figure [5l the time interval between panels a and b is 500 s, while that between 
panels b, c, and d is only 100 s. Thus, the evolution of the fiux-rope current structure 
exhibits a very distinct transition during this early development (see also the movie in the 
electronic version). Prior to time t' = 3.16 x 10^ s, the fiux rope rises in the corona only 
very slowly. At time t' = 3.17 x 10^ s and later, the fiux rope has strongly accelerated and 
is rising very fast. We evaluated the radial location and fiow speed of the plasma at the 
location of the 0-point during this interval, shown in Figure [61 Between t' = 3.07 x 10^ s 
(initial appearance of the 0-point) and t' = 3.16 x 10^ s, the fiux rope does not move 
significantly, remaining roughly at hpji = 1.2Rq with a speed of VpR — 20 km s~^. Between 
t' = 3.16 X 10^ s and t' = 3.24 x 10^ s, the fiux rope rapidly ascends in the corona as its 
velocity increases to Vfr{3.2Rq) ~ 900 km s"^ in 800 s. The average fiux-rope acceleration 
during this interval is apR = 1.1 km s^^. Thereafter, the CME moves out at a roughly 
constant speed of VpR ~ 900 km s~^ beyond hpR = 3.2Rq. Our simulation clearly yields a 
strongly accelerated, fast eruption. 

The evolution of our eru pting fiux rope is consistent with the three distinct phases 



described for observed CMEs (jZhang et al.ll200ll ): the initiation phase, when the fiux rope 



slowly rises in the corona; the rapid acceleration of the CME in the corona, corresponding to 
the impulsive phase; and finally the propagation of the CME in the interplanetary medium 
at a constant speed, the terminal phase. 

The fast eruption induces the formation of a new current sheet at the front of the fiux 
rope (panel c in Figure E]). This current sheet results from the formation of a compression 
region between the front of the erupting fiux rope and the ambient medium traveling at 
the speed of the isothermal solar wind. At the time of eruption onset, the plasma speed 
at the leading edge of the fiux rope is VLsi— 1.767?©) ~ 150 km s~^, significantly faster 
than the local solar wind speed, Vsw — 100 km s~^. Thus, the magnetic field is expected 
to be squeezed in this region, consistent with the development of a current sheet. However, 
comparing the fiux-rope velocity with the local Alfven speed, ca — 2000 km s~^, we conclude 
that this compression region is not a shock. 
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6. Interaction between the CME flux rope and the interplanetary medium 

The right column of Figure [5] shows the dynamics of the flux rope that propagates 
out of the corona and into the interplanetary medium after the eruption. The field lines 
representing the open flux of the northern polar coronal hole in Figure |5] are colored as 
a function of their component, with green and blue corresponding to \B^\ = G and 
\Brp\ = 0.02 G, respectively. The initial magnetic field (Figure |3]) is strictly poloidal, lying 
in the (r, 6) plane and having i?<^ = 0. The photospheric footpoint motions introduce a 
nonzero B^ into the system only within the central sheared arcade, below the south lobe of 
the fan. Part of this toroidal flux resides on the field lines that reconnect at the flare current 
sheet and become entrained into the disconnected flux rope If the ejected flux rope 

subsequently reconnects with an open line of the background field, that initially green line 
will acquire a nonzero B^ component and become partially blue. 

In panel d of Figure the innermost open green field line becomes blue all along a 
section that has a helical shape. The blue/green color and the morphology show that this 
previously open line has reconnected with the flux rope. The newly reconnected fleld line 
belongs both to the flux rope and to the coronal hole, coupling the flux rope to both the base 
of the corona within the polar hole and the remote heliosphere at the far end of the open 
field line. As the flux rope propagates through the northern coronal hole, this reconnection 
continues to process open flux, forming additional new fleld lines that connect the flux rope 
to the solar pole and to the interplanetary medium (panels e and f in Figure [5]). This 
third episode of reconnection, occurring between the flux rope and the open fleld, is similar 
to the interchange reconnection in a standard closed null-point topology. Since it involves 
closed CME fleld lines and open coronal-hole flux, we refer to it as the CME-interchange 
reconnection mode. 

Meanwhile, a fourth reconnection episode occurs. The fleld lines closed below the fan's 
north lobe (pink lines) and the open fleld to the south of the outer spine line (red lines) 
reconnect together. This reconnection is of the classical null-point interchange type, and the 
two resulting reconnected fluxes are the closed red arcade lines and the open pink flux shown 
in panels e and f of Figure [51 This episode is a direct consequence of the eruption of the flux 
rope and is due to the re-formation of the low-coronal null point (see ^|7]for details). After the 
red and pink fluxes have been processed fully in this way, reconnection across the re-formed 
null point continues between the yellow and orange flux systems. The latter reconnection 
is of the classical closed/closed type, and completes the restoration of the helmet streamer 
conflguration (Figure [5]F) from which the event began (Figure Hj). 

At the end of the simulation, the reconnection between the flux rope and the open fleld 
terminates well before all of the magnetic flux in the coronal hole has reconnected. The 
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magnetic flux of tlie flux rope is much smaller {^fr — 1-96 x 10^° Mx) than the open flux 
{^open — 2.26 X 10^^ Mx). Therefore, the flux rope entirely reconnects and disappears as a 
separate entity, while a great deal of unsheared open flux remains in the northern coronal 
hole. The total merger of a flux rope with the open interplanetary field through reconnection 
may occur during the eruption of faint CMEs with small magnetic flux. However, we do not 
believe that it should be the conventional evolution for all CMEs. For our simulation, we 
expect that the twist added to the open field lines will be dispersed through torsional Alfven 
waves well before the magnetic structure reaches the Earth. This would be in disagreement 
with the detection at Earth of interplanetary coronal mass ejec tions/magnetic clouds that are 
the c ounterpart of the CMEs launched during solar eruptions (jwimmer-Schweingruber et al. 
2006h . 

7. Evolution of the magnetic topology 

The reconnection dynamics described above give rise to a sequence of topological tran- 
sitions in the magnetic field. We note that throughout this paper, the words "fan" and 
"spine" have been used to denote the different separatrix structures of the magnetic field. 
It is important to keep in mind that these terms actually are defined for a 3D null-point 
topology, rather that the 2.5D null-line geometry of our simulation. In 3D, the fan lines form 
a surface while the spines are singular lines; in 2.5D, on the other hand, both structures are 
surfaces. Thus, there is no real physical distinction between the "fans" and "spines" in the 
simulation presented in this paper. We still employ those terms, however, in order to develop 
insight into the generalization of our results to fully 3D configurations. 

Figure [7] displays the six topologies that successively develop during the evolution. As 
before, separatrix field lines are plotted in dark blue, while other field lines are color-coded as 
in Figures m and |5] to show the same flux domains. Light blue field lines are added to Figure [7] 
to represent the sheared flux of the inner arcade and the twisted CME flux rope. The initial 
magnetic field corresponds to a closed null-point topology, associated with the null point 
NPl and its closed outer spine line (Figure [7^). The classical null-point separatrices delimit 
the different connectivity domains (see §2.3p . As described in ^ the subsequent null-point 
reconnection at NPl eventually opens the outer spine to the heliosphere, transforming its 
topology into the open null-point type (Figure [Tb). 

Following the onset of the eruption, the CME flux rope forms through flare reconnection, 
implying the pairwise formation of an 0-point 01 and a new X-type null point NP2 (Figure 
[Tt). 01 is located at the central axis of the flux rope, and NP2 is created in the flare current 
sheet below the flux rope, where the magnetic field lines of the sheared arcade are anti- 



- 16 - 



parallel. Accordingly, new separatrix surfaces associated with 01 and NP2 appear in the 
system and define new connectivity domains. The fiux rope belongs now to an independent 
fiux domain, delimited by a new separatrix SI in Figure [7t. The post fiare-reconnection 
loops below the erupting rope are delimited from the highly stressed original arcade by the 
separatrix S2, which is joined to SI by the current sheet containing NP2 (Figure [Tt). 

Although the two X-type null points, NPl and NP2, share some properties, we em- 
phasize that their physical origin and topological roles are very different. NPl is a robust 
topological feature due to the multipolarity of the field, and must be present independent 
of whether the system is ideal or dissipative. The breakout current sheet forms as a con- 
sequence of the deformation of NPl. NP2, on the other hand, is not a robust topological 
feature. First, the fiare current sheet forms as a result of the extreme stretching of the inner 
arcade field lines; then, NP2 and 01, and almost certainly many other null points, form 
inside this current sheet as a result of dissipation. NP2 and 01 would not occur under a 
truly ideal evolution and, in principle, they could disappear by simply merging with each 
other. In contrast, NPl can never disappear, even in 3D. Note also that the fiare current 
sheet associated with NP2, just like the breakout current sheet associated with NPl, has 
much more structure than is illustrated in the figures. During the highly dynamic phase of 
the simulation, both current sheets break up into multiple magnetic islands, each of which 
defines its own separatrix. These islands, however, eventually merge with the fiux at either 
end of the current sheets. Consequently, we show in Figure [7| only the time-averaged, global 
separatrix structure. 

During the early rising phase of the fiux rope (©, the current sheet due to the de- 
formed null point NPl becomes increasingly elongated, which increases the rate of breakout 
reconnection. Note that this reconnection acts to decrease the fiux inside SI, i.e., in the 
erupting fiux rope. On the other hand, the rate of fiare reconnection at NP2 increases even 
more rapidly, which results in a net increase in the fiux inside SI. This competition between 
the breakout and fiare reconnection continues until the two separatrices that can be seen at 
the bottom center of Figure [7t - the inner spine surface on the left and the fan surface on 
the right - converge and merge at the fiare current sheet NP2. 

The magnetic topology at this instant of merger is shown in Figure [TJi. Note that the 
topology is degenerate, in that all of the separatrices and the two null points are connected. 
This instant also corresponds to the disconnection of the CME fiux rope from the fiare loop 
arcade, at least in 2.5D. Consequently, the fiux in the ejecta stops growing at this time; from 
now on, that fiux only decreases, via CME- interchange reconnection at NPl. Note also that 
the null points NPl and NP2 switch nature at this time. NP2 now becomes the robust null 
point associated with the multipolarity, while NPl becomes the transitory null point that 
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can disappear by merging with 01 (Figure [Tjl). In fact, this is exactly what occurs later in 
the simulation. 

The start of the fourth reconnection episode disrupts the merged separatrices of 
Figure [Tli and forms a new topology. We use the label E for the new separatrices. The 
exchange of connectivity between pink and red field lines, at the null point NP2, transfers 
the open outer spine from the left (north) side of the flux rope to the right (south) side 
(Figure E^). Thus, the open null-point topology previously associated with the original null 
point NPl (Figure [7)d) now belongs to the flare-sheet null point NP2 and its inherited inner 
and outer spines (Figure [7^). The original null point NPl now lies on the separatrix surface 
El separating the disconnected CME flux rope from the background open flux. NPl acquires 
a spine that splits from the fan below NPl and the open outer spine above NPl (Figure [7^). 

The CME-interchange reconnection episode (^, which couples the flux rope to the 
open magnetic field, does not change the new topology. This reconnection occurs at the 
null point NPl and transfers open magnetic flux through the SI spine and separatrices to 
the far (south) side of the flux rope (Figure [7f). As indicated in the figure, although the 
reconnection does not alter the topology, it does transfer sheared light blue field lines from 
the disconnected flux rope to the open heliospheric magnetic field. 

8. Discussion 

To understand how fiare-accelerated particles can directly access the open interplanetary 
magnetic field during a solar eruption producing a CME, we investigated the dynamics of 
the magnetic reconnection when a flux rope forms, erupts, and propagates into the solar 
wind. We performed a 2.5D MHD simulation in spherical coordinates of an axisymmetric 
null-point topology, embedded in a helmet streamer formed by an isothermal solar wind. 

We have analyzed the detailed topological changes associated with a fast breakout CME 
erupting into an interplanetary magnetic field opened by the solar wind. The initial magnetic 
field is a null-point topology with a closed outer spine that is confined below the helmet 
streamer. Prior to the initiation of the CME, the magnetic energy builds up due to slow 
footpoint shearing motions imposed within the arcade below the null point. The excess 
energy so introduced inflates the sheared arcade, distorting the null point above it into an 
extended current sheet and eventually inducing magnetic reconnection there. After all of the 
closed flux below the streamer top and above the null point has reconnected, the initially 
closed outer spine opens into the heliosphere: the null-point topology now is surrounded by 
open field in the northern coronal hole. Additional null-point reconnection occurs gradually. 
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until the inner sheared arcade becomes highly distended and narrows at high altitudes to 
develop an internal (flare) current sheet. The onset of reconnection across this sheet causes 
the system to lose its equilibrium, rapidly forming a flux rope that erupts into the corona 
at high speed and propagates away into the interplanetary medium. During this evolution, 
additional episodes of reconnection occur that reform the open-spine null-point topology at 
low altitudes and, most important, couple the closed fleld lines within the CME flux rope to 
the open fleld lines of the interplanetary medium. 

The fast nature of our CME is noteworthy as an exception to the results of previous 
2.5D simulations of breakout CME eruptions into the solar wind fivan der Hoist et al.l 120071 : 
Zuccarello et al.ll2008l |2009| . l2012l : ISoenen et al.ll2009l ). The major difference between those 
studies and ours is the strength of the assumed background and active-region magnetic fields: 
we assumed fleld strengths that are larger by about a factor of 10, with a concomitant dif- 
ference by a factor of 100 in the plasma beta. Thus, it is not surprising that our strongly 
magnetically dominated corona yields a much more explosive, higher -speed eruption. The 
slow streamer-blowout simulation performed by IZuccarello et al.l ( 12012| ) nevertheless exhibits 
magnetic-reconnection dynamics similar to those discussed in this paper. After the eruption 
of their CME, which results from the breakout and flare reconnection, the flux rope recon- 
nects with the c l osed m agnetic fleld below the helmet streamer. The evolution described in 
Zuccarello et al.l ( 120121 ) follows the same topological changes as ours described in §71 except 
that their flux- rope separatrix surface SI remains closed within the helmet streamer, rather 
than reconn ecting out into the wind-opened corona. The results obtained with 3D MHD 
simulations (Ivan der Hoist et al.l l2009l : ICohen et al.l l2010t iLugaz et al.l l201ll ) should display 
an analogous sequence of changes as the CME propagates into the interplanetary medium. 
However, the evolution of the magnetic topology has not been clearly established in those 
studies. The complexity of the three-dimensional topological objects makes the determina- 
tion of the associated changes a far more challenging task than in the much simpler 2.5D 
geometry assumed in this paper. 

In the standard CSHKP model of eruptive flares, the flare-accelerated particles are 
trapped in the CME and do not have access to the open interplanetary medium (see §1]). 
However, this standard model does not consider the CME's interaction with and propagation 
into the ambient magnetic fleld. Our study demonstrates that the dynamics of magnetic 
reconnection, by inducing a speciflc topological evolution, can couple the magnetic fleld of 
a CME to the open interplanetary magnetic fleld. The resulting reconnected fleld lines 
are rooted at the solar surface within the pre-existing coronal hole, connect to the remote 
heliosphere, and pass through the CME flux rope, providing a path for the prompt escape 
of flare-accelerated particles to the Earth. 
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Energetic particles in the corona have a small Larmor radius and follow the reconnected 
magnetic field lines. Electromagnetic emissions (UV, X-ray, 7-ray) occurring when particles 
impact the denser chromospheric layer are localized at the footpoints of the reconnected 



field lines (see e. g. iDemoulin et al.lll997l : iMasson et al.ll2009al : iReid et al.ll2012[ ). Therefore 



it is well established that the reconnecting magnetic field effectively channels the particles 
energized by fiares. In the magnetic configuration of our MHD simulation, the particles 
initially trapped in the fiux rope field lines can access the interplanetary medium through 
the newly reconnected field lines coupling the open and closed fields. Therefore, energetic 
particles traveling initially along the CME fiux-rope field lines may be injected onto the 
newly reconnected field lines that open to the interplanetary medium. 

Energetic particles in the corona generate radio emissions. The coronal localization of 
these radio sources, combined with complementary observations, provides obse rvational diag- 
nosti cs on the loca tion of energetic e l ectron s and their related solar phenomena (jPick &: Vilmer 
20081 ). Recently, iDemoulin et al.l (120121 ) performed a detailed analysis combining radio, 
white-light, hard X-ray, and EUV observations of the CME on 2001 April 15. In this study, 
they highlighted specific radio emissions that correspond to the fiare reconnection occurring 
below the fiux rope. They also found that the CME fiux rope observed in w hite light was 



simult aneously observed at radio way e length s; radio CMEs we r e disco vered by lBastian et al. 



( 120011 ) and confirmed by iMaia et al.l ( 120071 ) . IDemoulin et al.l ( 120121 ) argued that the radio 
emission of the CME fiux rope is produced by energetic electrons, initially accelerated by 
the fiare reconnection and injected along reconnected field lines wrapping around the fiux 
rope. Finally, they showed that several radio sources appear during the post-eruptive phase, 
on the edges of the radio CME. They argued that these radio sources correspond to ener- 
getic electrons accelerated during magnetic-reconnection processes between the CME and 
the ambient magnetic field. In addition, these multiple episodes were temporally consistent 
with the injection of energetic ele ctron beams into the interplanetary medium derived from 
type-Ill radio-burst observations (jPick fc Vilmerll2008l ). 



Our numerical results on the reconnection dynamics of a CME propagating into the 
interplanetary medium are consistent with the observational results and interpretation of 
Demoulin et al.l (120 12[ ) concerning both the fiare reconnection and the development of 
multiple reconnection sites between the CME and the ambient field It is noteworthy 

that new injections of energetic electrons into the interplanetary medium occurred almost 
simultaneously with the inferred multiple reconnection episodes. This evidence strongly 
supports the new model for particle injection established in this paper, in which the coupling 
between the closed CME fiux-rope magnetic field and the open field of the solar wind provides 
a path for energetic particles to escape. 
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Our model magnetically reconfigures the corona through multiple reconnection episodes, 
providing theoretical support for a process a lr eady proposed to explain observations during 



solar energetic particle events. iKlein et al.l ( 120111 ) showed that western CME-less flares 



accelerating particles in the corona, do not produce any SEP events at the Earth. Thus, 
CMEs may be an essential ingredient for particle escape into the interplanetary medium. 
Those authors argued that, among other mechanisms, magnetic reconnection between the 
closed CME and the open interplanetary magnetic field may be required to inject particles. 
During CME-less flares, this coupling cannot occur, leading to particle confinement in the 
corona and no SEPs. 

During some SEP events, energetic particles can be detected farther than 90° from 



the fl aring active region f lLinlll970l : iKallenrode et al.lll992l : iNitta et al.ll2006l : iRouillard et al. 



201ll ). Usually, these in-situ measurements are ascribed to CME-driven shock acceleration 
that also injects energetic particles over a wide longitudinal range. The scenario described in 
this paper provides a new perspective on these observations. Indeed, the coupling between 
open and closed fields does not occur near the active region in our simulation, but far away 
from it. Rather than injecting energetic particles through a CME-driven shock over a wide 
longitudinal range, multiple reconnections between the CME (where flare-accelerated parti- 
cles are trapped) and open magnetic field can allow the particles to escape from locations 
remote from the original fiare site. This would be esp ecially true if the open field had a cor- 
ridor topology as described in lAntiochos et al.l ( 1201 ll ) , which ca n map to a large longitudinal 



arc in the hehosphere as suggested by several studies of SEPs (jPick et al.ll2006l : iKlein et al. 
20081 : bresing et al.l [2012! ). 



The primary conclusion from our study is that, in the presence of nearby open fiux, 
interchange-type reconnection between an erupting fiux rope and the open fiux naturally 
leads to the escape of fiare accelerated particles. Our results, however, raise major theoretical 
and observational issues. On the theoretical side, the most important issue is the extension 
of the model to 3D. A key shortcoming of our 2.5D simulation is that the erupting rope 
disappears completely as a result of reconnection with the open fiux. Figure [Sj whereas, 
ICMEs/magnetic clouds are observed in association with impulsive SEPs at 1 AU. However, 
in 3D the breakout and the interchange reconnection will occur at the current sheet formed 
from a true null point and is likely to show strong variation along the axis of the erupting fiux 
rope, which itself will be fully 3D. We conjecture that in this case interchange reconnection 
will not follow the one-to-one ordered reconnection imposed by the 2D geometry, but may well 
occur with some fiux from deeper inside the fiux rope, forming a complex ejecta entangling 
open and closed fiux. This conjecture will be tested with future fully 3D simulations. 



On the observational side, the most important issue raised by our model is the quantita- 
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tive distribution of flux required for flare-particle escape. Note that in order for the erupting 
flux rope to undergo interchange reconnection with open field, the overlying closed flux has 
to be removed. In the pre-eruption configuration shown in Figure [31 the amount of overlying 
closed flux is small compared to the amount of erupting flux; in other words, the coronal 
hole is "near" the eruption and, hence, interchange reconnection easily occurs. On the other 
hand, if the amount of overlying closed flux is large then interchange reconnection becomes 
highly unlikely. Therefore, comparison of our model with observed events and prediction of 
whether flare particles will or will not escape requires detailed knowledge of the amount of 
flux in the various coronal systems prior to eruption. In principle, such information can be 
obtained from high-resolution magnetograms, but again, fully 3D eruption models will be 
required, because concepts such as "overlying" and "near" become much more subtle in 3D. 
The 3D simulations, therefore, will provide the definitive observational tests of our model 
for flare particle escape. 
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A. Analytical expression for the dipole ring 

In cylindrical coordinates (^,0,2;), the vector potential for a toroidal ring of magnetic 
dipoles is A = z)e^, where is the unit vector along the axis. It is defined by the 

integral 



fJ'z^ - [fJ'za + fi^{z - b)] cos(x - 



z) = Bod'R^ I . ' Z . : : 3/2 ^x (Ai) 



[^2 + a2 + (^-6)2-2a^cos(x 



where BQ^fi^e^ + fiz^z) is the magnetic moment, with Bq its field magnitude and fi^, Hz its 
direction cosines respectively along the ^-axis and z-axis, such that + /i^^ = 1; a, 6, 
and X ci^^^G respectively the ^, z, and cf) positions of the dipole in the cylindrical coordinate 
system; and d and Rq are respectively the depth of the dipole below the surface and the 
solar radius. Introducing the new variable ip = (x — + 7r)/2, where x ~ represents the 
angular separation between the local contribution to the ring dipole and the location of the 
point where the vector potential is computed, we have 
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cos(x — 0) = cos(2'?/' — tt) = — cos(2'?/') = — 1 + 2 sin^ -i^. 



(A2) 



With this variable change, the vector potential becomes 



'0 



'^^^ [f^z{^ + cl) + f^^iz — b)] — 2 [fi^a + fi^{z — b)] sin^ ip 



;^2„,,l 3/2 



[(^ + ay + - - M sin^ V^] 



d?/^, (A3) 



after exploiting the symmetry of the integrand of Eq. (1A3P to reduce the integration range 
from [0,7r] to [0,7r/2]. Introducing the elliptic parameter k, 



4ae 



{^ + ay + {z-by 



(A4) 



we rearranged Eq. ( ]A3l) and reduced it (IGradshteyn fc RyzhikI Il965l ) to complete elliptic 
integrals of the first and second kinds, K{k) and E{k), 

'0 \/l — k sin^ ip 

'.7r/2 



dip\ 1 — k sin 



(A5) 



The following final expression for the vector potential of the dipole ring was obtained: 

U + af + {z-bff^[{i-af + {z-bY] 
{U + a)^i. + {z- b)^l^] E{k) - 2 [a/i, + {z- b)^l^] {K{k) + k-' [E{k) - K{k)]]^ (A6) 

We evaluated the poten tial numerically using numerical fits to the elliptic integrals given by 
Milne-Thompsoru ( 1l972l ). and then derived the magnetic field from B(r) = V x A(r). By 
superposing the solar background dipole field and the magnetic field created by the dipole 
ring (Eq. IA6p . for the parameter values specified in §2.3[ we obtained the 2.5D axisymmetric 
quadrupolar magnetic configuration shown in Figure It is straightforward to show from 
the above expression that, for the special case 6 = 0, a = i?© — d, and d <^ i?©, |B| = Bq 
at the solar equator (z = 0, = -R©), whether the dipole moment is tangent to the surface 
{^z = 1) or normal to the surface (//^ = !)• This result motivated the particular form chosen 
for the normalization constant in Equation lAll 
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Fig. 1. — The standard model for eruptive flares and its implications for particle trapping. 
The left panel shows the rising CME plasmoid represented by the green arch, along with the 
flare reconnection below this arch. The right panel shows that all the field lines resulting 
from the reconnection remain closed and, hence, do not allow any flare accelerated particles 
to escape. 



- 29 - 




Fig. 2. — Post-relaxation atmospheric conditions at = 27° for r G [l,20]i?Q. Top panel: 
radial solar wind speed; bottom panel: plasma beta. 
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Fig. 3. — Pre- and post-relaxation magnetic fields. Left and middle panels: the magnetic 
configuration before turning on the solar wind at t = s and after the system reaches 
its quasi-steady state at t = 2.1 x 10^ s, respectively. The radial magnetic field at the 
photospheric surface is color-shaded, with Bj. G [—10,20]. Gray lines display magnetic field 
that initially is totally closed (left panel), but subsequently opens at the solar poles (middle 
panel); dark blue lines display the magnetic separatrices. Right panel: a zoomed- in view 
of the null-point topology. Important topological features are labeled, and the location 
and direction of the photospheric forcing are indicated by black arrows on the photospheric 
boundary. 
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Fig. 4. — Reconnection at the null point. The temporal evolution of the </> component of the 
current density (left column) and selected magnetic field lines (right column) are shown after 
onset of the photospheric shearing motions. The radial magnetic field at the photospheric 
surface is color-shaded as in Figure [31 The cf) component of the current density is gray- 
shaded; black and white correspond to RqcP ■ V x B = —7 and +7, respectively. For the 
color scheme used to draw the field lines, see text. 
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Fig. 5. — Flare reconnection and flux-rope eruption. The radial magnetic field at the pho- 
tospheric surface is color-shaded as in Figures [3] and HI and the component of the current 
density is gray-shaded as in Figure HI Field lines are drawn as in Figure HI except that the 
northern coronal-hole field lines now are colored according to the value of their component 
of the magnetic field: green and blue correspond, respectively, to \B^\ = G and \B^\ = 0.02 
G. (An animation is available in the online journal.) 
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Fig. 6. — Flux-rope trajectory. The top panel displays the temporal evolution of the height 
of the 0-point of the flux rope; the bottom panel displays its speed. 
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Fig. 7. — Schematic of the topology evolution during the CME eruption and coronal prop- 
agation. Light blue lines represent inner-arcade field lines with shear; the other field lines 
follow the same color scheme as in Figures S] and [51 Panels illustrate: a) the initial null-point 
topology with closed outer spine; b) the null-point-reconnected topology with the outer spine 
opened; c) formation of the fiux rope and its related separatrices due to fiare reconnection; 
d) merging of the separatrices and the formation of a transient complex topology; e) splitting 
of the merged separatrices to form two new open-spine null-point topologies; and f) transfer 
of magnetic shear from the closed fiux rope to the open interplanetary field. 



